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Calcium sulfate devolatization during the impact at Chicxulub, Mexico and dispersal in the strato- 
sphere of the resultant sulfuric acid aerosol have been suggested as a possible mechanism for the 
Cretaceous- Tertiary extinctions [1, 2]. In this paper, we investigated two shock-induced devolatiza- 
tion reactions of calcium sulfate up to 42 GPa in the laboratory: 

CaSC>4 + Si02 — + CaSiC>3 + SO3 "f 

CaSOf * CaO *4- SO2 j + j O2 T 

We found both to proceed to a much less extent than calculated by equilibrium thermodynamic calcu- 
lations. Reaction products are found to be 10 -2 times those calculated for equilibrium. Consequently 
our estimate of the amount of sulfur oxides degassed into the atmosphere from shock devolatization of 
CaSC>4 in the Chicxulub lithographic section (6xl0 15 -2xl0 16 g in sulfur mass) is lower by a factor of 
70 to 400 than previous estimates, the related environmental stress arising from the resultant global 
cooling of **4 K and fallout of acid rain does not appear to suffice to explain the widespread K-T 
extinctions. 
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1. Introduction 

Since Alvarez et a/.[3] proposed the Cretaceous-Tertiary (K-T) impact hypothesis, the Chicxulub impact 
structure has been recognized as the leading candidate for the crater(s) associated with the bolide. Strong 
supporting evidences relating the Chicxulub, Yucatan crater with the world-wide extinction include the 
following: 

I* Tektites, presumably representing fused target material, from Beloc, Haiti and Mimbral, Mexico K- 
T sections show the glass source terrane is likely to be a continental margin [4, 5]. This environment 
contains sediments rich in carbonate and sulfate, which is consistent with Chicxulub stratigraphy. 
Oxygen and strontium isotope analyses [6] also favor Chicxulub over the Manson, Iowa crater as the 
source for the tektite-Iike glassy ejecta, which, in most world- wide locations, has been weathered 
to various clays; 

2. The tektites and Chicxulub melt rock have been ( 39 Ar/ 40 Ar)-dated and found to be coeval within 
0.1 m.y. at 65 m.y.[7j; 

3. Tsunami deposits exactly at K-T boundary, which have now been recognized to occur around the 
present the Gulf of Mexico. This supports the hypothesis of a nearby impact site[8, 9]. 

The dominance of calcium carbonate (calcite) and calcium sulfate minerals (anhydrite and gypsum) in 
the upper 3 km of the Chicxulub section (in the sequence of limestone, anhydrite-dolomite conglomerates, 
limestone-dolomite and anhydrite-gypsum layers [10]) has led to active research into the effect of their 
devolatization upon impact. Previous shock recovery experiments conducted on single crystal calcite 
found calcite devolatization upon release to ambient pressure from shock pressures greater than ~ 10 
GPa. These experiments have suggested that a Chicxulub-size crater may be capable of increasing the 
CO2 budget of the atmosphere by 5xl0 18 to 2 xl 0 19 g, therefore the global atmospheric inventory of CO2 
(2.2xl0 18 g in pre-industrial age [2], 2.7xl0 18 g at present [11]) could have been increased by a factor of 
2 to 10. This is hypothesized to have caused global warming of 2 to 10 K for periods of 10 4 to 10 5 years 
[11], An environmental stress, which affected the proliferation of the wide range of genera which became 
extinct at the K-T boundary has been suggested to have been induced by impact-liberated SO2 or SO3 
from sulfates. Estimates of the amount of SO2 and SO3 devolatization (based on equilibrium conditions) 
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ranges from 4x10 to 8xI0 18 g. If this amount of sulfur oxides produced stratosphere H2SO4 aerosol, 
a rapid temperature decline of 10 to 20 K is predicted to have occurred for several years [1, 2]. In this 

paper, we present results of calcium sulfate devolatization experiments under shock conditions via two 
reactions: 

CaS04 + SiC>2 — +■ CaSiOa 4- SO3 f 

CaSC>4 — * CaO + SO2 | + 5 O2 t (B) 

Reaction (A) is of interest because of its lower incipient reaction temperature than reaction (B) (1200 K 
vs. 1450 K) based on Gibbs formation energy calculations using thermodynamic data in [12] (we assume 
ambent SO2/SO3 and 0 2 partial pressures of 10~ 4 and 0.2 bar respectively as in the contemporary 
atmosphere). Heating experiments carried out in platinum tubes at 1 atm and 1, 047-1, 406°C suggest 
this materia] could be the origin of the yellow Haiti glass[4]. The basement of Chicxulub is generally 
thought to provide quartz for the reaction, but it is not clear at present how anhydrite/gypsum at 
shallower depths could have reacted with the basement rocks. Nonetheless, the two reactions have been 
investigated with anhydrite/silica powder mixture and anhydrite/gypsum polycrystalline disks. 

Comparing to theoretical calculations, we found the experimental degree of sulfate devolatization for 
powdered samples (particle size was tens of microns) was - 10 2 times lower than theoretically estimated. 
Also, we estimate the mass of sulfur oxides released during the formation of the Chicxulub crater using 
the upper limit of experimental reaction efficiency. Moreover, upon calculating the anisotropic shock 
pressure decay in the lithologic rock section in which the Chicxulub crater formed, we estimate is that 
6xl0 15 to 2xl0 16 g of sulfur in various oxides could have been generated. Finally, we discuss the possible 
environmental effects of injecting this amount of sulfur in the form of sulfuric acid and oxides into the 
stratosphere. 

2. Recovery experiments 

Solid recovery experiments were performed on the Caltech 20 mm and 40 mm propellant guns. Samples 
(~ 6 mm in diameter and ~2 — 3 mm in thickness) were encased in stainless steel containers, which were 
pumped to ~ 30 millitorr vacuum before the shots [13, 14]. 
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2.1. Reaction between silica and anhydrite 

Experimental parameters are listed in Table 1 (two shock experiments on anhydrite/SiO are listed and will 
be discussed later). The starting material was a mixture of silica (crystalline or amorphous, Alfa #13024 
and 89709) and natural crystalline anhydrite (Ward’s Geology #46E0535) powders and was pressed into 
target container to initial densities of 60-85% of its Archimedian density. Average silica grain size was 4 
fim, and anhydrite grains were mostly between 30-100 ^m. Equation-of-state constants for the mixture 
(see Table 2) were calculated from previous anhydrite data of Simakov et a/. [15] and quartz data of Swegle 

ei a/.[16] using the formulae by Boslough [17]: assuming uniform stress distribution, for a two-component 
system, 

V = rmVi + 7712^2, (1) 

Kos = [(vi/Kosi) + (t>2/A’os2)]“ 1 > ( 2 ) 

KqS — -^0s[ v l(l "h •^■0Sl)/-^0Sl d - v 2(l + -^ 052 )/^ 0 S 2 ] — 1) (3) 

where mi, Vi, Vi, Aosi and A '' 0Si are the components’ mass fractions, specific volumes, initial volume 
fractions, bulk moduli and their pressure derivatives at zero pressure. Shock pressures determined by 
impedance-match method range from 27.4 to 42.3 GPa. Five 20 mm shots and the two 40 mm shots 
(see Table 1) all had 1:1 molar ratio of anhydrite : silica. The initial porosity of the mixture were 
17.2% (1106) to 40.0% (1108); Shot 1107 employed fused quartz; The 40 mm shots were conducted to 
determine the shock duration effect on the reaction, but no reaction was seen in shot 917, and shot 923 
was not recovered; Finally, shots 1109 and 1110 had anhydrite : silica molar ratios of 1:3.7 and 3.0:1. 
Recovered samples were analyzed with petrographic microscopy, scanning electron microscopy (SEM, 
instrument: Camscan Series 2 with TVacor Northern EDS detector TH-3/54-6901, operated at 15 kV) 
and X-ray diffraction (XRD, instrument: Scintag DMC-008, radiation source: Cu-K Ql ). Compared with 
the original material, the changes exhibited in the 20 mm (except shot 1111, which will be described 
separately in the following paragraph) post-shock samples are quite similar: in agreement with previous 
research [18], silica becomes amorphous in spite of its original crystallinity (see XRD spectra in Figure 2); 
Anhydrite is recovered as a crystalline phase. Although shock-induced mosaicism in the crystal grains 
was observed with cross-polarized light on petrographic microscope, it appears unlikely that anhydrite 



5 


recrystallized from a melt as no rounding of the grains was observed (see Figure 1 a). 

In shot 1111 (in which 10 mass% iron powder was intentionally mixed in addition to anhydrite/quartz), 
devolatization was much more extensive than the rest of the shots and reaction of iron to iron sulfate 
and iron sulfide were observed. In their study of sulfur speciation in basaltic glasses [19], Carroll and 
Rutherford reported that proportion of dissolved sulfur present as sulfate (as opposed to sulfide) increases 
from near 0% at FMQ (fayalite-magnetite-quartz) oxygen fugacity to near 100% at 2 to 3 log /o 2 units 
above FMQ. The oxygen fugacity present in these recovery experiments was well into the sulfate stable 
regime. The presence of iron sulfide led us to believe the greater degree devolatization of CaSO-j within 1 
mm of the stainless steel container in all the 20 mm shots was affected by the reducing effect of the metal 
and would not have occurred in its absence. In the central metal-free region, the dimensions of possible 
reaction zones are so limited that they were nearly at the limit of spatial resolution of the SEM. In the 
following three sections we will attempt to derive the actual degree of devolatization from experimentally 
observed chemical compositions at different locations in the samples. 

2.1.1. SEM instrument resolution 

The SEM electron beam spot is much less than l^m, but the dimension of excitation volume in the sample, 
and therefore the instrument resolution, is larger due to electron scattering and secondary fluorescence 
in the sample [20]. A “smearing” function is assumed to convolve with the “true” chemical composition 
to give the observed composition. The function form is taken to be Gaussian: 


/(*) 



(4) 


/ f(x)dx = 1 (5) 

OO 

where 2 1 is a measure of spatial resolution. Anhydrite and quartz disks of ~ 1 mm thickness each were 
sandwiched together and heated at 573 K for 6 hours, followed by a 24-hour press at 4000 psi so that 
plastic flow may take place to produce a good contact (with less-than-1 //m gap) as a no-reaction reference. 
A comparison between the SEM analysis across the pressed boundary and calculation (convolution of 
Equation 4 and a step function) found / = 0.53/im to provide the best fit (Figure 3). This agrees with 
our expectation that the resolution distance is larger than the electron beam diameter. 
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The data in Figures 3 and 4 are corrected for secondary fluorescence excited by characteristic ra- 
diations. Another concern has been the fluorescence excited by the continuous spectrum. We use an 
approximate equation ((15.10) in Reed [21]), modified for compounds by mutiplying the ratio of mass 
attenuation coefficients of the excited element A (/. t *) and the compound {nc), the intensity of fluo- 
renscence Ij relative to electron-excited characteristic K-radiation of element A in the compound I £ is 
thus 


4- = 9-7x10 - 6 Z 4 ^, 
l c He 


( 6 ) 


(Z is the atomic number of the excited element). The correction factor for continuum fluorescence is 

Fj = 1/(1 + 4 ). (7) 

L c 


We consider three cases: 


1. Compound is 50 mol.% CaS 04 and 50 mol.% Si02 (in the middle of the mixing zone), 
Ij a /Ig a = 1-21%, If/1%= 0-24%; 

2. Compound is near 100% Si02 with a trace of Ca and S (deep into silica). a / l^ a ■=. 0.67%, 

7//j£=0.13%; 

3. A sharp CaS 04 -Si 02 boundary, with the electron beam shifted to Si02 side so that the electron- 
excitation volume is completely in Si02- Any Ca and S signal is purely due to secondary fluorescence 
from the Si02 continuum (Si characteristic line is not energetic enough to excite Ca or S). The 
secondary fluorescence intensity, relative to electron-excited radiation in pure CaS0 4 , is given by: 

T Ca < S Cfl,5 

+S- = 0.5 x 9.7 X 10-e^'^Oi. (8) 

*CaS0 4 HCaS0 4 

Equation 8 is very similar to Equation 6 except the factor 0.5, which arises because only half of 
the continuous radiation goes into CaSC> 4 , neglecting the finite width of the primary X-ray source. 
The ratios for Ca and S calculated are 0.23% and 0.11%. 

More detailed numerical calculations were done and the results agree within ±0.3%. Case (3) agrees very 
well with observations at 1 //m in quartz from the cold-pressed boundary, where Ca and S signal intensities 
are 0.2±0.1% and 0.4±1% of those in pure CaS0 4 (Figure 3). For the shock-recovered samples, the Ca 
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and S intensities within ^ 3/rni from the boundaries are much higher than the secondary fluorescence 
level and the corrections are negaligible compared to the analytical uncertainty. 

2.1.2. Mixing of sulfate with silica 

In shocked samples the boundary layer between calcium sulfate and silica is thicker than the cold-pressed 
edge. In the following we examine several possible mixing mechanisms: 

1. Solid state diffusion: 

The diffusion constants of H, 18 0 and 30 Si atoms in quartz have been documented in [22]. At 800°C, 
they vary over a wide range, with H having the highest D = 2.5 x 10" 11 m 2 /s, and 30 Si having 
the lowest D = 1.3 x 10“ 21 m 2 /s. In the time scale of our experiments ~ lfis, the characteristic 
distance y/Dt ~ 10‘ 8 -10~Vn is much smaller than the observed reaction zone thickness; 

2. Liquid state diffusion: 

Ruble et al. directly measured oxygen self-diffusivity in Na 2 Si 4 0 9 melt up to 1825°C and between 
4-20 GPa [23]. The diffusion constant they reported ranges from 1.0 to 4.2xlO" 10 m 2 /s increasing 
with temperature and pressure. Si-0 bond breaking is the basic process controlling both O self- 
diffusion in Na 2 SL|0 9 and CaS0 4 diffusion in silica melt which is of present interest. It is possible 
to use these data to obtain an order-of-magnitude estimate of mixing time and length scales in 
the Si0 2 liquid. Again, in the Ifis shock duration, the highest diffusivity (4.20xl0 _1 ° m 2 /s) gives 
characteristic distance 2 x 10- 2 /zm, which is still too small to account for the /mi-size mixing layer; 

3. Rayleigh-Taylor instability: 

Although initially crystalline quartz is amorphized during shock, we cannot conclude that it has 
been once molten because quartz can transform to diaplectic glass without melting [18]. Since there 
is a strong contrast in strength of quartz and anhydrite (~ 1 GPa for quartz and ~ 0.1 GPa for 
anhydrite), we suggest Rayleigh-Taylor instability as the third mixing mechanism. 

Rayleigh-Taylor instability arises at interfaces between two materials of different strength when they 
are strongly accelerated or decelerated along a direction perpendicular to their planar interface. 
According to the theoretical model by Drucker[24], when shock wave propagates from the stronger 
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material into the weaker material, the interface experiences alternating compressional and tensile 
stress due to perturbations (bumps) on material surface. When the stress difference exceeds the 
strength of the stronger material (<jq), the bumps grow freely and instability occurs. Two important 
derivations of the theory are threshold perturbation amplitude: 

W = H{ 1 + */2)v 0 F/P, (9) 

and time dependence of instability growth above the threshold: 

h-hl h = (h 0 - hl h ) cosh y/0Pt/(\pH) t (10) 

where H is material thickness, ho is initial perturbation amplitude, p is material density, P is 
shock pressure, A is the perturbation wavelength, F is a geometric factor between 0.25-1, and 

/* = rriV~2-s. 

A prominent feature of the theory is that the threshold is independent of wavelength A. Experiments 
by Barnes et ai on aluminum and 304 stainless steel plates support the theoretical prediction [25]. 

For our case, o’o — 1 GPa, p=2.65 g/cm 3 , P=40 GPa, A — H ~ grain size 4pm, the threshold 
thickness perturbation calculated from Equation 9 is ^ 0.4pm, which is very reasonable. Above 
threshold, the growth is very fast: for an initial perturbation 10% above threshold, Equation 10 
indicates that it takes a few’ nanoseconds to grow to sizes comparable to A/2, after which the theory 
is no longer valid. 


From the above discussion, only Rayleigh-Taylor instability emerges as a plausible mixing mechanism. 
Although it is not a diffusive process, we describe the shocked sample interface mathematically by a linear 
diffusion profile. Denote gca(x) “the molar ratio of Ca/O, the solution of one dimensional diffusion 
equation (heat equation, see, e.g., [26]) with initial conditions: 


9Ca(x)\ t= o - 


l 


0.25 

0 


if x < 0 
if x > 0 


is 


9Ca(x) = gerfc(x/I). 


(U) 
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A fit to the shocked sample profile (convolved with /(*)) gives the mixing length L ~ 1.5/tm (solid 
line in Figure 3): 

9ca M ( x )=J f{x - x')g C a{x')dx' (12) 

2.1.3. Degree of devolatization 


The experimentally seen devolatization shows some scatter (Figure 4). For simplicity, we assume an 
average of 20% S molar loss in the region -X < * < X, so the ratio S/O is 


0s(z) — O.Sgca(x) = 0.1erfc(x/X). ( 13 ) 

The S/Ca profile observed on SEM is given by 

/ S/r w x _ J- 00 9s(x , )e-( t '-* ) 3 f l7 dx' 

gca(x')e-i‘'-*)’H'dx'' ( 14 ) 

Bulk devolatization, defined as the fraction of sulfur loss in the boundary layer over the total sulfur 
mass in the original sample, is given by 


DV = 3 X f—L^yCa ~ 9s)dx 

9cAr—^R ’ ( 15 > 

where R is the anhydrite grain size. The factor 3 takes into account the three dimensional effect.For 

R = 100pm, evaluation of the formula yields a numerical value of DV = 6 x 10 -3 . 

We infer that anhydrite mixed into silica upon shock loading and underwent devolatization during 

release. Tyburczy and Ahrens [13] used the following approach to calculate the extent of shock-induced 
reactions: 


1. Entropy excess required for incipient reaction: 



where T 0 is room temperature, T,r is the temperature of incipient reaction (at which the sums of 
Gibbs formation energies for the reactants and products are equal), C r is the atmospheric pressure 
heat capacity at constant pressure. 

2. Entropy excess required for complete reaction: 

Scr = Sir + AS- ^ riiRln(Pi/P 0 ) l (17) 

gas products 
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where AS is the entropy difference between reactants and products as computed from Robie et al 
[12], the last term on the right takes into account effects of partial pressures, P,, of gas products 
(Po is the ambient pressure), n, is the number of moles of gas specie t, R is the gas constant (8.31 
J mol'^K" 1 ). 

3. Entropy gain in the shocked state (and in the post-shock state assuming isentropic release): 

ASh = S tr + C v \n(Tjj/Ts), (18) 

where S ir is the entropy change of phase transition during compression, 7s is the temperature 
of isentropic compression from initial volume Vo at temperature 7} to Hugoniot volume V H and 
isentropic pressure P s : 



where 7 is the Griineisen parameter. 7# is the shock temperature, determined by: 

Vh f TH 

-(P H -P,) = J C v dT ( 20 ) 

C v is heat capacity at constant volume (at high pressure). 

4. The extent of reaction is given by: 

Fraction of material reacted = — ~ — (91 'i 

Scr - Sjr ^ 

The results from Equation 21 for reactions (A) and (B) are shown as the solid and broken curves in 
Figure 5. The equation-of-state parameters are the same as in Table 2, and St r =0 for the calculation. 
Local devolatization in the reaction layer is close to theoretical calculation, but decomposed anhydrite 
is only a small portion (~ 6 xl 0 “ 3 ) of the total mass because of the dimension of the reaction layer is 
much smaller relative to the grain size (~ 100 /zm). In the Chicxulub section, mixing of calcium sulfate 
and quartz was certainly less intimate, the reaction efficiency should be even lower. 

2.2. Decomposiiion of anhydrite and gypsum alone 

Since calcium sulfate strata often do not occur close to quartz in the Chicxulub section, decomposition of 
anhydrite/gypsum is an important issue to the evaluation of sulfur oxides production. Similar recovery 
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experiments up to 42.7 GPa have been conducted on anhydrite and gypsum samples (Table 3). The 
samples were nonporous single crystalline gypsum and polycrystalline anhydrite disks approximately 4 
mm in diameter and 0.5-1 mm in thickness (mass 20-30 mg). In several of the shots the samples were 
sandwiched between Zn disks in an attempt to enhance solid recovery. 

On SEM images (Figure 1 b), the recovered samples show textures that are not different from un- 
shocked anhydrite and gypsum. There is no indication of inhomogeneous shear band heating which has 
been recognized as the major volatile loss mechanism for calcite [27]. Thermogravimetric analysis (TGA) 
was performed on Thermal Analysis System SETAEAM TG92. Milligram- mass samples were analyzed in 
a He atmosphere, heating from 25 to 1600°C at a heating rate of 2°C/min. Impact-induced mass loss was 
determined as the difference in mass loss between shocked and unshocked samples. TGA and chemical 
analyses on the SEM (on apparatus described previously) revealed no variation in sulfur concentration 
beyond analytical uncertainty (~ 2%). Devolatization of 72-76% of the total H 2 0 was detected by TGA 
in shocked gypsum (Table 3). 

In summary, in the pressure range covered by the experiments, the fraction of sulfur, by mass, that 
degassed is ~ 10“ 2 of theoretical prediction (Figure 5), and overall sulfur devolatization is minimal. 
Both reactions (A) and (B) are endothermic, with rather large enthalpy increases of 313 (for o-quartz 
as reactant, 305 for fused quartz) and 502 kJ per mole reacted CaS0 4 [12], which are 17 and 27% 
of the specific shock energy of a 40 GPa impact on non-porous materials. To examine the role of 
enthalpy difference in shock-induced reactions, we conducted two recovery experiments on anhydrite and 
amorphous silicon monoxide (SiO, Alfa #89430) powder mixture (Table 1). The reaction in question is 

CaS0 4 + SiO — * CaSi0 3 + S0 2 T (C) 

We observed extensive reaction in the recovered material of shot 1098, where S in anhydrite is reduced 
by a factor of 4-5. Although enthalpy of SiO glass is unknown, it is a less stable compound than quartz 
and therefore absorbs less energy in reaction (C). The weaker bonding structure compared to Si0 2 could 
also give rise to a much higher diffusion rate. Both may contribute to the excessive reaction occurred. 
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3. Implications for the K-T boundary 

S.l. Pressure decay in the transient Chicxulub crater 

Previous studies on crater formation mechanisms give formulae for pressure decay with distance R along 
the centerline in a transient crater[29, 30, 31]. For example, Ahrens and 0’Keefe[29] calculated gabbroic 
anorthosite or iron meteorites cratering planetary anorthosite target and fit their calculations with two 
branches using exponential functions of the form: 


kgioWMbar) = a log 10 (R/R 0 ) + b (22) 

where i?o is the radius of the meteorite, parameters a and b assume different values for different targets 
and in the near- and far-field (a < 0). 

The Schmidt-Holsapple scaling [31] breaks down the pressure decay in large craters into four ranges- 
near source, strong shock, intermediate shock and material strength regimes, in each of which the peak 
shock pressure decays exponentially with normalized distance R/R 0l the exponents (equivalent to pa- 
rameter a in Ahrens-0 ’Keefe model) are 0, -3.6, -1.8 and -1.18 in successive regimes away from impact 
center. 

These two approaches give roughly the same pressure decay profiles if the parameters are chosen 
properly. A more important effect to take into account in evaluating sulfate devolatization is the de- 
viation of equi-pressure contours from a hemispherical surface because of rarefaction wave from the 
target free surface. Define Victual an< * ^ actual to be vaporization and melt volumes (normalized to 
impactor volume) given in O’Keefe and Ahrens [30], and V£f lt to be the normalized vaporiza- 

tion and melt volumes calculated assuming hemispherical peak shock pressure contours. For a 15 km/s 
anorthosite meteorite impacting anorthosite V r ” t e Ji;/V r A m / ,t is only 0.052; For a 45 km/s anorthosite on 
anorthosite impact, V^ t e t [ t o j/V; m / /t =0.040, V r a ^ r ,/V^“ por =0.042; For a 15 km/s iron on anorthosite im- 
pact, VZcfl'ai 0.426. In this paper, we modify the simple model by assigning each hemispherical 
shell an “average” pressure that is lower than the centerline pressure. Although this is not a rigorous 
treatment, it can be regarded as a first-order correction for the effect of the two-dimensional rarefaction 
wave (Figure 6). The pressures of the points labeled “complete vaporization” and “complete melting” are 
complete phase change pressures, their radii are chosen such that the hemispherical volumes within are 
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equal to the volumes given in [30]. Such derived a, 6 parameters for the “average” pressure attenuation 
are listed in Table 4. 

A few noticeable features in Figure 6: 

1. Near impact center, pressure profile is dominated by the initial shock, centerline and average pres- 
sures should not differ very much; 

2. Their deviation grows larger at larger radius as rarefaction wave becomes important; 

3. At great depth (far away from the target surface), self-similarity is established, the average pressure 
curve is parallel but offset downward to the centerline pressure curve. 

In Figure 7, centerline pressure profiles given by Holsapple-Schmidt scaling [31] and Ahrens-O’Keefe 
formula [29] are plotted together with estimated averaged pressure for transient Chicxulub crater, assumed 
to behave as gabbroic anorthosite, if the bolide was an chondritic asteroid, also similar in behavior as 
gabbroic anorthosite. The effect of different initial impact velocity is represented through the center 
peak pressure. The average pressure decay is then deduced from the Ahrens-O’Keefe profile and the 
corrections in Table 4. Based on this profile, the theoretical incipient devolatization pressure (~ 25 GPa) 
corresponds to a devolatization hemispheric diameter of 50 km in the target rocks. 

An alternative scenario is the bolide was a comet (p ~ 1 g/cm 3 , t; ~ 60 km/s). From in siiu 
mass spectrum measurements in several Halley flybys by the Giotto and Vega spacecraft in 1986, it 
was found [32] the element composition of the comet is similar (within a factor of two) relative to Cl 
chondrites for elements heavier than O. The comet is much more enriched in CHON elements, and as a 
result, the atomic concentrations of heavy elements are about one fourth of the corresponding values for 
chondrites. Since the average atomic weight is also lower for the comet (~ 9.6 vs. ~ 14.9 for chondrites), 
the comet/chondrite mass abundance ratio for heavy elements is about 0.4. To maintain the global Ir 
anomaly constant (estimated at 3.7xl0 10 g from the Gubbio section [3]), the comet’s diameter has to be 
around 14 km. Following a similar calculation as for a gabbroic anorthosite asteroid, the devolatization 
hemispheric diameter would become 97 km. 
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3.2. Conclusions 

Our experiments do not cover the entire pressure range of the K-T impact, but our results show conclu- 
sively that shock-induced calcium sulfate devolatization requires higher shock pressure than devolatization 
of serpentine, calcite and Murchison carbonaceous chondrite[13, 27, 33]. The extent of devolatization re- 
actions are on the order of 10“ 2 less than predicted via equilibrium calculation. One reason for this 
appears to be the lack of shear bands in anhydrite/gypsum and poor mixing between anhydrite and 
silica. Scaling up to the Chicxulub crater gives a rough estimate of sulfur degassing during its formation. 
Assume: 

1. A 500 m-thick calcium sulfate bed 1 km below surface [1] in which 50% of the mass was CaS04, 


2. the degree of reaction was about 2% (upper limit of experimental uncertainty of anhydrite/gypsum 
devolatization) within the reaction diameter, 

the total degassed sulfur mass wais 6xl0 15 g (in forms of SO2 and SO3) for an asteroid impact. Under 
the same assumptions, the estimate in case of a comet impact was 2xl0 16 g. 

This mass is much less than previous estimates given by Brett[l] or Sigurdsson et al. [2] (4xl0 17 - 
8xl0 18 g) mainly due to the experimentally observed devolatization occurs to a much lower extent (by a 
factor of 70 to 400) than previously assumed, but it is still about one to two orders of magnitude above 
the amount of SO2 released by the eruption of Tambora in 1815 (estimated at 1-2 xlO 14 g and caused 
cooling of 0.7 K in the following 2 years [34, 36, 37]). However, our estimate supports Blum et al.’s 
suggestion that most of the Haiti yellow glass was derived from shock-induced mixing of silicates with 
carbonates rather than sulfate-bearing rocks[6]. 

How the reactions proceed at pressures under which the reactants are molten is not answered by our 
present experiments. To give an upper limit, we assume perfect mixing and 100% devolatization inside 
34 km hemispheric diameter from the center of the Chicxulub crater (which corresponds to pressures 
greater than 52 GPa), the amount of S released is about 1.4 xlO 17 g, which is still a factor of 3 lower 
than previous estimates. However, we regard such an amount as unlikely. 
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3.3. Environmental effects 

The theory of lofting of impact volatized SO2 and SO3 from a crater and temperature decline due to the 
increase in optical depth resulting from stratospheric H2SO4 aerosol has yet to be worked out. Therefore, 
evaluation of the environment effects of the K-T impact mostly relies on climatic changes after historic 
large volcanic eruptions. Sigurdsson used the following empirical relation to correlate the volcanic sulfur 
mass yield x, in grams, to observed temperature decreases after the eruptions [35] 


AT = 5.9 x 10~ 5 x 0 ' 31 (23) 

Extrapolating to our estimate of sulfur release (t=6x 10 15 -2x 10 16 ), it is expected to give rise to 
short term global cooling of -4-6 K and acid rainout of 40-150 g/m 2 on the earth’s surface. Recently, 
Vogelmann et al. [38] discussed the effect of volcanic SO2 release on the ozone layer by providing H2SO4 
aerosol surfaces on which chlorine can react to forms that catalyze ozone depletion [39, 40]. Vogelmann et 
al. point out that the resulting increase in UV-B radiation (290-320 nm) could be biologically hazardous 
even for large volcanic eruptions. However, H2SO4 aerosol also increases optical depth in the UV radiation, 
which could partially compensate for ozone depletion. The relation between S0 2 release and ozone 
depletion is not established at present. In conclusion, we suggest that in a few years after the K-T impact, 
the volatized sulfur oxides and sulfuric acid aerosol would have caused a drop in the earth’s temperature 
by 4-6 K, but the cooling could have been partially or completely offset by the C0 2 green house effect. 
As sulfuric acid aerosol separated from the atmosphere, the earth underwent a more persistent warm 
period of tens of thousands of years when the temperature was raised by as much as 10 K [11]. The 
environmental effects of sulfate devolatization could have been quite hazardous, but itself alone may not 
suffice as a potent extinction mechanism. 
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Tables 


TABLE 1. Recovery experiments on anhydrite/silica and anhydrite/SiO. 


Shot 


Sample 


Molar ratio (anhydrite:silica) Projectile (km/s) P pe>k (GPa) 


1106 

1107 

1108 

1109 

1110 
1111 
1112 

917(40mm gun) 
923(40*nm gun) 

1098 

1099 


anhydrite/qtz 
anhydrite/fused qtz 
anhydrite/qtz 
anhydrite/qtz 
anhydrite/qtz 
anhy drite/qtz/Fe 
anhydrite/qtz 
anhydrite/qtz 
anhydrite/qtz 
anhydrite/SiO 
anhydrite/SiO 


1:1 

1:1 

1:1 

1:3.7 

3.0:1 

1:1 

1:1 

1:1 

1:1 


Ta 

W 

Ta 

Ta 

Ta 

Ta 

Ta 

Ta 

Ta 

W 

Ta 


P r °j 

1.87 
2.02 
1.92 

1.90(?) 

1.94 

1.89 

1.88 
1.59 
1.77 
2.06 
1.58 


33.8 

42.2 

35.0 

42.3 
32.5 

34.3 

34.1 

27.4 

31.5 


1 On Caltech’s 20mm gun unless noted otherwise. 



TABLE 2. Equation-of-state constants of 
anhydrite, silica and their mixture. 


Material 

P (g/cnr 3 ) 

A'os (GPa) 

KU 

Anhydrite (LPP) 

2.97 

38.5 

6.0 

Silica (HPP) 

4.29 

350 

3.3 

Mixture (1:1 molar) 

3.28 

48.7 

7.6 
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TABLE 3. Recovery experiments on anhydrite and gypsum. 


>hot 

Sample 

Projectile 

Pproj (km/s) 

P pe»k (GPa)^ 

H 2 O loss (%) 

1090 

anhydrite 

Ta 

1.11 

23.6 



1091 

anhydrite 

Ta 

1.48 

33.2 


1092 

anhydrite 

Ta 

1.66 

38.0 

- 

1103 

anhydrite 

Ta 

1.77 

41.2 



1093 

anhydrite 

Ta 

1.82 

42.7 



1101 

gypsum 

Ta 

1.35 

29.7 

72 

1094 

gypsum 

Ta 

1.59 

36.2 

76 

1102 

gypsum 

Ta 

1.81 

42.1 

76 


1 On Caltech's 20mm gun. 

2 Calculated using stainless steel equation of state [28], 



TABLE 4. Centerline and average pressure attenuation: 
lo 8io(-P/Mbar) = a log 10 (.R/flo) -f 6 




Near-field 

Far-field 



a 

6 

a 

b 

An-* An 
45 km/s 

Centerline 1 

-0.293 

1.208 

-2.15 

2.373 

Average 

-0.616 

0.947 

-2.18 

1.394 

An— >An 
15 km/s 

Centerline 1 

-0.222 

0.285 

-1.97 

1.399 

Average 

-0.608 

-0.024 

-1.97 

0.556 

Fe— *An 

Centerline 1 

-0.295 

0.470 

-2.49 

1.475 


15 km/s Average -2.49 1.167 

1 After Ahrens and O’Keefe [29]. 



TABLE 5. Parameter s a, b used to estimate peak and average pressure s of Chicxulub crater. 

Near-field Far-field 

a b a 5 

Centerline -0.222 0,501 -1.97 1.62 

Average -0.608 0.192 -1.97 0.772 
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Figures 



Fig. 1. SEM images of after-shock samples: a. Anhydrite (light contrast grains) and quartz (dark area in between). Shot 
#1106, 33.8 GPa; b. Anhydrite/gypsum. Shot #1102, 42.1 GPa. 
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Fig. 2. X-ray diffraction spectra of silica-anhydrite mixture before- and after-shock (Shot #1109, 42.3 GPa). The two 
spectra at bottom are JCPDS standards for CaSO* and SiC> 2 . Molar ratio of silica to anhydrite is 3:1. Initially crystalline 
quartz is amorphized in the after-shock material. 
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Fig. 3. Atomic ratio of Ca/O profiles of three recovery shots on calcium sulfate and quartz. A cold pressed sample was 
used to estimate SEM resolution. Theoretical fits are given by Equation 12 with parameters: (1) dashed curve: l = 0.53/im, 
L—0, (2) solid curve: l = 0.53 (im, L = 1.5Aim. The data at 1 ^m on the dashed curve represent three analyses at different 
locations near the boundary. 
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Fig. 4. Atomic ratio of S/Ca profiles of recovery shots on calcium sulfate and quartz. Solid curve is given by Equation 14. 
Further into Si0 2 , S/Ca ratio becomes indeterminate and is indicated by the dashed line. Error bars represent SEM 
analytical uncertainty. Results for other 20 mm shots are similar and not shown. 
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Fig. 5. Degree of devolatization via reactions (A) and (B) as a function of shock pressure. The curves are calculated using 
Equations 16-21. Partial pressures of 0 2 and S0 2 /S0 3 are taken to be 0.2 and 10" 4 bar respectively in the calculation, 
which are representative values for normal atmosphere. Impact velocities are indicated for shock pressures of 20 40 60 and 
80 GPa. 
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F»g. 6. Shock pressure versus normalized radius at 45 km/s for gabbroic anorthosite impactors. Holsapple-Schmidt scaling 
relation is also shown for the same impedance-matched center peak pressure. Inset: Pressure attenuation for: a. 15 km/s 
anorthosite on anorthosite impact; b. 15 km/s iron on anorthosite impact. See text for meaning of ’’average". 
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Fig. 7. Shock pressures versus normalized radius for transient Chicxulub crater from an asteroid impact. The bolide are 
target rock are assumed to have equation of state of gabbroic anorthosite. Bolide initial velocity is 20 km/s, Ro—5 km. 
Parameters used are listed in Table 5. See text for meaning of "average”. 
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